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. Abstract. The PHENIX experiment at RHIC has measured transverse energy and charged 

particle multiphcity at mid-rapidity in Au + Au collisions at y^s^^iv ~ 19.6, 130, 62.4 and 200 
GeV as a function of centrality. The presented results are compared to measurements from 
other RHIC experiments, and experiments at lower energies. The y^sTriv dependence of dEr/drj 
and dNch/drj per pair of participants is consistent with logarithmic scaling for the most central 
^ , events. The centrality dependence of dEr/dri and dNch/drj is similar at all measured incident 

fT^ ' energies. At RHIC energies the ratio of transverse energy per charged particle was found 

CNj I independent of centrality and growing slowly with y^s^^iv"- ^ survey of comparisons between the 

. data and available theoretical models is also presented. 
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1. Introduction 

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Q ■ Laboratory was designed to measure the properties of matter at extremely high temperatures 

^ , and densities. Under such conditions, there is possibility to produce states of matter that have 

yet to be observed or studied in the laboratory. The best known of these is the quark-gluon 
plasma (QGP), a form of matter where quarks are not confined within individual baryons but 
^ I exist as some form of plasma of individual quarks and gluons. It should be emphasized that the 

' exact properties of this matter are not known and that the characterization of the deconfined 

state, if such a state is produced, is an essential part of the RHIC program. 

One fundamental element of the study of ultrarelativistic collisions is the characterization 
of the interaction in terms of variables such as the energy produced transverse to the beam 
direction and the number of charged particles. These two variables are closely related to the 
collision geometry and are important in understanding global properties of the system during 
the collision. 

This paper describes a systematic study of dEx/dr] and dNch/drj at mid-rapidity^ by the 
PHENIX experiment at center-of-mass energies ^/s^ = 19.6, 130, 62.4 and 200 GeV. The 
centrality dependence of Et and Nch is characterized by the number of participants, determined 
with a Glauber model, and is studied as function of the incident energy. Et and Nch results for 
all four RHIC measurements are included as part of this study. The data taken at 19.6 GeV is 
particularly interesting because it allows a close comparison to the lower energies of the CERN 
SPS program. Comparisons are also made to previous experiments at the Brookhaven AGS and 



^ For the full list of authors and acknowledgments see reference Q 
^ Also referred as Et and Nch in this paper. 



CERN SPS at center-of-mass energies of 4.8 GeV, 8.7 GeV, and 17.2 GeV. Finally, an extensive 
set of collision models describing the Et and Nch distributions are compared to the existing 
data. 

2. PHENIX detector and Analysis 

PHENIX is one of four experiments located at RHIC The PHENIX detector consists of two 
central spectrometer arms, designated east and west for their location relative to the interaction 
region, and two muon spectrometers, similarly called north and south. Each central spectrometer 
arm covers a rapidity range of |?7| < 0.35 and subtends 90*^ in azimuth. The muon spectrometers 
both have full azimuthal coverage with a rapidity ranges of —2.2 < i] < —1.2 (south) and 
1.2 < r] < 2.4 (north). Additional global detectors are used as inputs to the trigger and for 
global event characterization such as vertex, time of event and centrality determination. A 
detailed description of the PHENIX detector can be found in [S]. 

The PHENIX detector subsystems relevant for the physics analysis published here are: the 
Pad Chambers used for the charged particle multiplicity measurement, the Electromagnetic 
Calorimeter used to measure transverse energy, the Beam-Beam Counter and the Zero Degree 
Calorimeter, the two detectors used for triggering and centrality determination are described in 
publications H El El 13. 

The analysis procedures to measure Et and N^h are described in details in publications piU 
1111 1121 [T^ . Some additional information concerning the analysis at ^/s^ = 19.6 GeV can be 
found in pQ. The preliminary results on dNch/dr] at ^/s^ = 62.4 GeV use the same technique 
as implemented at other energies. 

For the sake of space we only mention errors relevant to the analyses. Statistical errors are 
small and do not exceed 1% of the measured value. The systematic errors are summarized in 
Tabled The errors for the lowest and the highest incident energy are listed, whereas at all other 
energies the are close to these values. The systematic errors for both measured values are of 
two types. The first type affects the centrality dependence. It is listed in the Tabled with the 
range (first number corresponds to the most central bin). Errors of the second type contribute 
to the overall scaling of the data. In the figures below the centrality dependent errors are shown 
as a corridor, and the the full systematic error is shown with the error bars^. Centrality related 
errors are common to both measured quantities. 



Table 1. Summary of systematic errors given in %. When the range is given, the first number 
corresponds to the most central bin and the second to the most peripheral bin. 

dEr/dT] dN^h/dr] 



V^NN [GeV] 


19.6 


200 


19.6 200 


Energy resp. 


4.7 


3.9 




Bkg. / noise 


0.5-3.5 


0.2-6 


1.0 1.0 


Acceptance 


2.0 


2.0 


2.3 2.3 


In- &; outflow 


3.0 


3.0 


5.7 2.9 


Occupancy 






1.6-0.3 3.5-0.1 


Centrality 


2.0 


0.5 


same 


Np 


2.9-6.7 


2.8-15. 


same 


Trigger 


0.4-8.8 


0.3-16. 


same 



^ Here and everywhere errors correspond to one standard deviation. 



3. Results 

3.1. PHENIX results 

The distribution of the raw transverse energy, Et^mc •> fiducial aperture of two EMCal 

sectors is shown in the left three panels of Fig. ^ for three RHIC energies ^s^^ = 19.6, 130 
and 200 GeV. The lower scale corresponds to the fully corrected Et normalized to one unit of 
pseudorapidity and full azimuthal acceptance. The lower axis in the plot is not labeled beyond 
200 GeV to avoid confusion between the true shape of the dEx/dr] distribution and Et as 
measured using the limited acceptance of two EMCal sectors. 




dET/dTi [GeV] ° 200 dN,,/dTi 



Figure 1. The distribution of the raw Et in two EMCal sectors (left) and the number of tracks 
in the east arm of the PHENIX detector (right) per MB trigger, measured at three energies. The 
lower axis corresponds to mid-rapidity values of dET/drj and dNch/dr] respectively. Distributions 
of the four 5% most central bins are also shown in each plot. 

Two EMCal sectors each with azimuthal coverage = 27° were used to make a measurement 
at ^Sj^j^ = 130 GeV. At other energies we used five EMCal sectors. Results obtained with 
difi^erent numbers of sectors at the same energy are consistent within 1.5%. 

The right three panels in Fig. ^ show the number of tracks reconstructed in the east arm 
of the PHENIX detector after the background subtraction and all other corrections. The lower 
axis corresponds to measured distributions normalized to one unit of pseudorapidity and full 
azimuthal acceptance. For a similar reason as for the Et , the lower axis is not labeled above 
200 GeV in dN^h/dq. 

For the Nch measurements at ^/s^ = 130 GeV, only the east arm was used, while for the 
other energies the measurements were made using both PHENIX central arms. The results 
obtained with two arms at ^/s'^ = 200, 62.4, and 19.6 GeV are consistent with each other 
within 1.5%. 

The distributions shown in Fig. ^ have a characteristic shape with a sharp peak that 
corresponds to the most peripheral events. Missing events caused by the finite MB trigger 
efficiency in peripheral events would make this peak even sharper than measured. The plateau 
in all distributions corresponds to mid-central events and the fall-off to the most central Au + Au 



events. The shape of the curves in Fig. ^ in the fall-off region is a product of the intrinsic 
fluctuations of the measured quantities and the limited acceptance of the detector. 

The distributions for the four most central bins 0%-5% to 15%-20% are also shown in each 
panel. The centroids of these distributions are used to calculate the centrality dependence of Et 
and Nch ^- The statistical uncertainty of the mean values (less than or about 1%) determined 
by the width of the distributions are small because of the large size of the event samples. 

The magnitude of dEx/dr] at mid-rapidity divided by the number of participant pairs as a 
function of Np is shown in Fig. |2 The right three panels show the same ratio for dNch/dr] at 
three RHIC energies. 
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Figure 2. dEx/dr] (left) and dNch/dij (right) divided by the number of participant pairs at 
different RHIC energies. Errors shown with vertical bars are full systematic errors. Lines show 
the part of the systematic error that allows bending or inclination of the points. Horizontal 
errors denote the uncertainty in determination of Np. PHENIX preliminary result for Nch at 
^/s^ = 62.4 GeV is shown in the right panel. 



The horizontal errors correspond to the uncertainty in Np, determined within the framework 
of the Monte Carlo Glauber model. The vertical bars show the full systematic errors of the 
measurements added quadratically to the errors of Np. The lines denote the corridor in which 
the points can be inclined or bent. The statistical errors are smaller than the size of the markers. 
The upper panel also shows the results of the two lower panels with open markers for comparison. 

An important result for Fig. El is the consistency evident in the behavior of the centrality 
curves of Et shown on the left and N^h shown on the right for all measured energies. Both values 
demonstrate an increase from peripheral (65%-70%) to the most central events by (50-70)% at 
RHIC energies ^/s^=130 GeV and 200 GeV. For the lowest RHIC energy (^/s^=19.6 GeV) 
this increase is at the level of systematic uncertainties of the measurement. One can note that 
results from PHOBOS jSIl) show that the total charged particle multiplicity is proportional 

* All plotted and quoted numbers correspond to the centrality bin-by-bin average values or ratios of corresponding 
averages. 



to Np while the multiphcity at mid-rapidity over Np increases with Np, indicating that the 
pseudorapidity distribution gets more narrow for central events. FigurelSJalso show the PHENIX 
preliminary result for Nch at ^/s^ = 62.4 GeV added to the central panel. 

The ratios of the dEx/dr] and dNch/drj per participant pair measured at different RHIC 
energies are shown in Fig. |21 In these ratios some common systematic errors cancel. 
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Figure 3. Ratios of dET/dr] (left) and dNch/dij (right) measured at different RHIC energies. 
The errors shown with vertical bars are the full systematic errors. Lines show that part of the 
systematic error that allows bending or inclination of the points. The horizontal errors denote 
the uncertainty in the determination of Np. 

The increase in the Et production between 19.6 GeV and 200 GeV (with an average factor 
of 2.3) is larger than for Nch (with average factor of 1.9). This is consistent with an increase in 
the particle production per participant common to both Et and Nch and a ~20% increase in 
(my) of produced particles contributing to the Et parameter only. 

The ratio of 200 GeV/19.6 GeV shows an increase from peripheral to central events; however 
the increase is marginal in comparison to the systematic errors of the measurement. 

The ratio of 200 GeV/130 GeV is flat above Np ~ 80 and is equal to 1.140 ± 0.043 for Et 
and 1.126 it 0.036 for Nch in the most central bin. A rather sharp increase in the ratios of 
both quantities between Np =22 and 83 in the ratios of both quantities is still at the level of 
systematic uncertainties. 

The preliminary ratio of 200 GeV/62.4 GeV although shows the decrease with Np which is 
within the systematic errors. The increase in the particle production between these two energies 
is about 40%. 

The ratio of the transverse energy and charged particle multiplicity at mid-rapidity as a 
function of centrality is shown in Fig. |l]for the three energies. The upper plot also shows the 
results displayed in the lower panels for comparison. 

The ratio ET/Nch ^, sometimes called the "Global Barometric Observable", triggered 



^ Er/Nch is used as a shortcut for {dEr / drf) / {dNch/ dr)) at = in C.M.S. 
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Figure 4. Er/Nch vs. Np at different RHIC 
energies. Tfie errors shown with vertical bars 
are the fuh systematic errors. Lines show 
that part of the systematic error that allows 
bending or inclination of the points. The 
horizontal errors denote the uncertainty in the 
determination of Np. 



considerable discussion |2H1 HH]- It is related to the (my) of the produced particles and is 
observed to be almost independent of centrality and incident energy of the collisions within the 
systematic errors of the previous measurements. The present paper forges a direct link between 
the highest SPS and lowest RHIC energies, making a more quantitative study of Et/Nch possible. 

The results presented in Fig. |3 show that the centrality dependence of ET/N^h is weak and 
lies within the systematic errors plotted with lines. There is a clear increase in Et/Nch between 
y^sj^=19.6 GeV and 200 GeV. The ^/s^ dependence of the results is discussed below. 



3.2. Bjorken Energy Density 

The Bjorken energy density [30] can be obtained using 



eBj 



1 dET 
A±T dy 



(1) 



where r is the formation time and A± is the nuclei transverse overlap area. 

The transverse overlap area of two colliding nuclei was estimated using a Monte Carlo Glauber 
model A±_ ~ (Jxf^y, where Gx and Uy are the widths of x and y position distributions of the 
participating nucleons in the transverse plane. The normalization to ttB? , where R is the sum 
of r„ and d parameters in a Woods-Saxon parameterization (see e.g.: |l]), was done for the most 
central collisions at the impact parameter 6 = 0. 

For the transformation from dET/dr]\ri=o to dEx / dy\y=Q, a scale factor of 1.25 it 0.05 was 
used, see P for the details of conversion procedure. 

The Bjorken energy density for three RHIC energies is plotted in the left panel of Fig. [21 For 
the 5% most central collisions, esj-T was 2.2ib0.2, 4.7±0.5 and 5.4ib0.6 GeV/(/m^ -c) for ^Js^ 
=19.6, 130 and 200 GeV, respectively. These values increase by 2%, 4% and 5%, respectively, 
for the maximal A'paj.j=394, as obtained from extrapolation of PHENIX data points. There is 
a factor of 2.6 increase between the "SPS"-like energy (y^s^ =19.6 GeV) and the top RHIC 
energy ^/s^ =200 GeV. The comparison of the only published eBj=3.2 GeV/fm^ at SPS |HJ 
and top RHIC energies, assuming the same r=l fm/c, reveals an increase in energy density by 



a factor of only 1.8, which may come from an overestimation in the SPS measurement, as shown 
in the left panel of Fig. Ilfll below. 

Another approach is used by STAR in for the estimate of the transverse overlap area 
of the two nuclei A± ~ Np^^ in Eq. ^ This approach accounts only for the common area of 
colliding nucleons, not nuclei. The results differ only in the peripheral bins as shown in the right 
panel of Fig. For a comparison, the same panel shows the result obtained by STAR which 
agrees with PHENIX result within systematic errors, displaying a smaller increase of the energy 
density with Np. 
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Figure 5. cbj • t deduced from the PHENIX data at three RHIC energies (left) and using 
different estimates of the nuclear transverse overlap at ^/s^ =130 GeV (right). 

Can we derive an estimate for r? We might say, on general quantum mechanical grounds, that 
a particle of full energy close to m-r can be considered to have "formed" after a time t = h/rriT 
since its creation in that frame. 

To estimate the average transverse mass, we use the final-state Et/Nch results presented in 
Fig. 0J As it was mentioned above, this value is basically independent of centrality and changes 
from 0.74 GeV to 0.88 GeV for ^Js^ from 19.6 to 200 GeV. Such estimate would provide us 
with the formation time r which is almost the same for all initial conditions of the collisions at 
RHIC. 

An approximate factor of 2/3 is used to account for the total number of particles making the 
estimate of (mr) ^0.57 GeV at ^/s^ =200 GeV and 130 GeV and {mr) =^0.50 GeV at ./s^;^ 
=19.5 GeV. This results in a formation time of r ~0.35 fm/c and 0.40 fm/c at these energies 
respectively. Then for the 5% most central events the estimate of the Bjorken energy at the 
three RHIC energies are 5.5 GeV//m^, 13 GeV/fm^ and 15 GeV / fm^ . 

3.3. Comparison to other measurements 

Several factors complicate the comparison between the results of PHENIX and the results of 
other experiments. AGS and SPS data were taken in the Laboratory (Lab.) system while the 
RHIC data are in the Center of Mass (C.M.S.) system. Since r/ and Et are not boost invariant 
quantities, the data should be converted into the same coordinate system. Some experiments 
provide a complete set of identified particle spectra from which information about Et and N^h 
can be deduced, while other experiments require additional assumptions to extract dEx/dr] and 
dNch/dr] from their results. Publication [J discusses it in details 

The PHENIX results for N^h are compared to the data available from the other RHIC 
experiments. This comparison is shown in the left panels of Fig. |21 
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Figure 6. Left panel: dNch/dij per pair of A'p measured by the four RHIC experiments at 
different energies. The shaded area is the PHOBOS systematic error. Right panel: RHIC 
average values (including PHENIX) compared to the PHENIX results. 



There is good agreement between the results of BRAHMS (221110], PHENIX, PHOBOS [SI 
I42[ I43j and STAR 44 , 45^ using Np based on a Monte-Carlo Glauber model. This agreement is 
very impressive because all four experiments use different apparatuses and techniques to measure 
the charged particle production. The systematic errors of all results are uncorrelated, except 
for errors due to the same Glauber model which are small. That makes it possible to calculate 
the RHIC average and reduce the systematic uncertainty. The averaged results from all four 
RHIC experiments are plotted in the right panel of Fig. El The procedure of the calculating the 
average is the same as used in [22 and is also explained in 

Figure [7| compares Et results from the PHENIX and STAR jlHl experiments. The results 
are consistent for all centralities within systematic errors, though STAR dEx/dr] per participant 
pair as shown in the upper panel has a smaller slope for Np in going from from semi-peripheral 
to central collisions and Ex/N^h shown in the lower panel is consistent for all Np. 

The RHIC run at y^s^=19.6 GeV allows us to make a connection between RHIC and SPS 
data. The highest SPS energy of 158A GeV corresponds to ^/s^ = 17.2 GeV in the C.M.S., 
making a direct comparison of RHIC and SPS results possible. This comparison is shown in 
Fig. H The SPS data is taken from [HI (111101 HnHMlilBl details of the data compilation are 
explained in pQ. 

Several comments should be made about this comparison. For both measured parameters 
the PHENIX results and the SPS results agree. The WA98 results are systematically higher, 
especially for dEj-fdr]. However the WA98 data has an additional systematic error common to 
all points shown for the last bin. For Net the relative spread of the SPS results is larger than 
for the RHIC results shown in Fig. (01 though overall the y^s^=17.2 GeV SPS measurements 
are consistent with the PHENIX result at ^s^^^ = 19.6 GeV. The NA57 results at this energy 
and at lower energy are published without systematic errors so they cannot be compared 
to other results on the same basis. 
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Figure 7. dEx/dr] divided by the number 
of Np pairs (top) and Et/Nch (bottom) 
measured by the PHENIX and STAR [MI 
experiments at ^/s^ =200 GeV. PHENIX 
systematic errors are explained in the text. 
The shaded area is the STAR systematic 
scaling error. 
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Figure 8. dEj-fdri (left) and dNch/dt] (right) divided by the number of Np pairs measured 
by PHENIX at y^s^^=19.6 GeV (solid markers) and recalculated from the results of the SPS 
experiments at the highest energy y^s^=17.2 GeV (open markers). The p + p result of NA49 
is marked with an open cross. 



Different SPS and AGS experiments made measurements at lower energies. The combined 
data of AGS, SPS and RHIC provide a complete picture of the centrality behavior of Et and 
Nch as a function of the nucleon-nucleon energy. The centrality dependence of dNch/drj at 
mid-rapidity measured at ^s^j^=4.8, 8.7 and 17.2 GeV by different experiments is shown in 
Fig. M The data is taken form [13 EOi EZl EHl EHl IZOl EH IZl ESI Ei IZ3 and the details of data 
compilation can be found in pp. 

At the highest SPS energy the averaging procedure is the same as for RHIC energies and 
weighted experimental errors are scaled with the reduced x^-like factor S reaching the value of 
1.5 at some points. For the intermediate SPS energy ^Sj^^=8.7 GeV, two experiments NA45 |37j 
and NA50 [3H1 reported the centrality dependence of dNch/drj at mid-rapidity. The discrepancy 
in the measurements is approximately three times the quadratic sum of their systematic errors. 
However the shapes of the two curves are almost the same. NA49 has published results |191 120j 
which give one point in dNch/drj at Np =352. This point favors the NA45 result. The NA57 
results also plotted in the Fig. |^ demonstrate even larger discrepancy with the NA50 
measurements, but these points lack the systematic errors, as we mentioned before. 

The average centrality curve is produced by taking into account the shape of the centrality 
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Figure 9. dNch/drj divided by the number 
of Np pairs measured by AGS and SPS ex- 
periments and the average taken at different 
energies recalculated in the C.M.S. 



curves reported by NA45 and NA50 and the single NA49 point. The errors are scaled with the 
factor S, which reaches a value of 2.5 at some points. The AGS results are presented with a 
curve produced from the combined results of the E802/E917 experiments (see publication ll ). 

The average SPS centrality dependence at ^/s^ = 17.2 GeV shown in the upper panel in 
Fig. ini and the average curve of the two RHIC experiments at ^s^^=19.6 GeV shown in the 
lower panel in Fig. El are very similar. Less than 5% increase is expected due to the difference in 
the incident energies between the highest SPS and the lowest RHIC energies (see section [3.4.11 
below). 

The average values presented in Figs. El and IH] are summarized in Table [2 



3.4- Dependence on the incident nucleon energy. 

The data compilation made in the previous section allows for a detailed study of the charged 
particle production in heavy ion reactions at different incident energies of colliding nuclei. 
Although the data on transverse energy production is not abundant, a similar comparison can 
be made [TOIIII]. 

3.4-1- Central Collisions Figure ITfll shows the energy dependence for the dEx/dr] and dNch/drj 
production per pair of participants in the most central collisions measured by different 
experiments. See ^ for the details of the data compilation. 

The results shown in Fig.^Jare consistent with logarithmic scaling as described in |in[ ll3 | [T2]. 
Use of the logarithmic function is phenomenological and is suggested by the trend of the data 
in the range of available measurements. The agreement of the fits with the data in both 
panels is very good, especially in the right panel where the averaged values are used for Np 
=350. The single point of NA49 [SJ is excluded from the Et fit. The results of the fit 

dX/d7] = {0.5Np ■ A)ln{^SNN/ \Js%n) are: 

for Et ^/s%^ = 2.35 ± 0.2 GeV and A = 0.73 ± 0.03 GeV 



Table 2. Average values of dNch/drj / (O.bNp) at different ^/s^. An additional 5% error should 
be added to columns 17.2 GeV through 4.8 GeV for the uncertainty related to recalculation to 
the C.M.S.. 
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Figure 10. Left panel: dErp/dr] divided by the number of Np pairs measured in the most central 
bin (value given in brackets) as a function of incident nucleon energy. The line is a logarithmic 
fit. The band corresponds to a \a statistical deviation of the fit parameters. Right panel: the 
same for dNch/dr]. The values of N^h are the average values corresponding to Np = 350. The 
single point at ^Sj^^ = 56 GeV is based on [IT] . 

for Nch \/4n^ = 1.48 ± 0.02 GeV and A = 0.74 ± 0.01. 

The parameter ^ s^j^=2.35 GeV obtained from the Et fit is slightly above, although within 
3(7 from the minimum possible value of ^/s^ = 2 x a. m.ii. =1.86 GeV. The measurement closest 
to it at ^/s^ = 2.05 GeV done by the FOPI experiment allows to estimate the amount of dEx/drj 
produced to be 5.0 GeV in the most central collisions corresponding to Np =359. Publication P 



gives details of the estimate. This does not disagree with the extrapolation of the fit but does 
indicate that the logarithmic parameterization requires higher order terms to describe how the 
Et production starts at very low ^/s^. 

The right panel of Fig. shows the logarithmic fit to the Nch data. It agrees well with 
all dNch/di] results plotted for A'p=350. Unlike that for Et, the fit parameter ^ for N^h 
is 1.48ib0.02 GeV which is lower than the minimum allowed ^/s^■ This suggests that above 
2 X a.m.u. the Nch production as a function of ^/s^ should undergo threshold-like behavior, 
unlike the Et production which must approach zero smoothly due to energy conservation. 

The FOPI measurement at y^s^=1.94 GeV and 2.05 GeV agrees with the extrapolation of 
the fit at energy very close to 2 x a.m.u.. It is an interesting result that colliding nuclei with 
kinetic energies of 0.037 GeV and 0.095 GeV per nucleon in the C.M.S. follow the same particle 
production trend as seen at AGS, SPS and RHIC energies. 

A fit to the charged particle multiplicity shows a factor of 2.2 increase in dNch/dij per 
participant in the most central events from the highest energy at the AGS i^/s^ =4.8 GeV) 
to the highest energy at the SPS {^/s^ =17.2 GeV) and a factor of 2.0 from the highest SPS 
energy to the highest RHIC energy (y's^ =200 GeV). Assuming the same behavior extends to 
the LHC highest energy y/s^=5500 GeV one would expect dNch/dr) = (6.1 ± 0.13) • (0.5A'p) 
and the increase in particle production from the highest RHIC energy to be ~60% for the most 
central events. With the greater energy, the rapidity width should increase by ~60% i.e. the 
total charged particle multiplicity would increase by a factor of ~2.6 from the top RHIC energy. 

It is interesting to compare the ^/sj^ dependence shown in Fig. to the similar plot for 
dNch/dy vs. ^/s^ published in In Fig. 4 of that publication the energy dependence of 

dNch/dy for Np =350 changes the slope between AGS and lower SPS energy. That is not seen 
in Fig. ^1 The parameterization suggested in publication |32], dNch/dy oc (^/stvjv)"'^ is also 
different from the logarithmic scaling in Fig. 1101 Although quantities in these two figures are 
not identical, some contribution to the shape of the curve in Fig. 4 in |^ may be due to use of 
logarithmic scale on Y-axis, producing a bend-like shape out of straight line in a linear scale. 

The ratio of Et /Nch for the most central bin as a function of -y/sj^ is shown in Fig. IllL Note 
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Figure 11. Ratio of Et over Nch for the 
most central events function of y^s^ 

recalculated into C.M.S.. The line is the 
ratio of two fits shown in Fig. 1101 The band 
corresponds to one standard deviation of the 
combined error. 



that the line shown in the figure is not the fit to the data points. Rather, it is calculated as a 
ratio of the fits plotted in Fig. IIOL The calculation agrees well with the data. 

There are two regions in the plot which can be clearly separated. The region from the lowest 
allowed ^/s^ to SPS energy is characterized by a steep increase of the ET/Nch ratio with ^/s^. 
In this region the increase in the incident energy causes an increase in the {rriT) of the produced 
particles. The second region starts at the SPS energies. In this region, the Et/Nch ratio is 
very weakly dependent on ^/s^. The incident energy is converted into particle production at 
mid-rapidity rather than into increasing the particle (w-t)- 



The shape of the Et/Nch curve in the first region is governed by the difference in the \J s^j^j^ 
parameter between Et and A^c/i- The second region is dominated by the ratio of the A parameters 
in the fits; this ratio is close to 1 GeV. Extrapolating to LHC energies one gets a ET/N^h value 
of (0.92±0.06) GeV. 

The Enp/Nf-h parameter is usually associated with the chemical freeze-out temperature T. A 
simple relation like 2/2>Et /Nch ~ mo + 3/2T finds T ~0.22 GeV for a measured value Et/N^.]^ of 
0.88 GeV at top energy at RHIC. Coefficient 2/3 takes into account neutral particle contribution 
and the average particle mass rriQ ~ 0.25 GeV can be estimated from jl7| . The freeze-out 
temperature estimate of T ~0.22 GeV extracted this way is in an agreement to T ~0.17 GeV 
measured from the spectra ^7] plus some contribution of particle flow. In publications OH 
authors suggest the same Et/Nck value close to 1 GeV is valid even at AGS and SIS energies. 
This is not in disagreement with the curve in Fig. ^2 because the Et definition used in |33| I34j , 
includes full baryon rest mass. We use different definition of Et as explained in ^ where mass 
of pre-existing baryons is not included. At lower energy where such contribution to the Et is 
larger than at RHIC it would bring the Et/Nch ratio closer to one. 

The transition between two region with different Et/Nch behavior occurs in the lower SPS 
energies. There are other observables which also undergo transition between regimes in the same 
range of ^/s^. For example, the {K^) / (tt^) ratio demonstrates a peak at around ^/s^ ^8 GeV 
as shown in Fig. 4 in reference |35j . In the same region, the volume of the colliding system 
goes through a minimum as shown in Fig.l in reference The asymmetric flow parameter 

V2 lies in-plane at and above the highest SPS energy and out-of-plane below it, see Fig. 8 in 
publication |2Z|- These evidences suggest that the energy range of 5-10 GeV has interesting 
physics potential which needs to be studied. 



3.4-2. Centrality shape Another interesting question is how the shapes of the centrality curves 
of Et and Nch change with ^/s^. 

One approach, previously used in a number of papers is to describe the shape of the centrality 
dependence as a sum of "soft" and "hard" contributions such that the "soft" component is 
proportional to Np and the "hard" component to the the number of binary collisions N^. 
A X Np + B X Nc. A disadvantage of this approach is that the contributions called "soft" and 
"hard" do not necessarily correspond to the physical processes associated with these notations. 
Another approach is to assume that the production of Et and N^h is proportional to N^; 
although the parameter a does not have any physical meaning. 

We present the results of B/A and a obtained from the fits to the data at different -^/sj^ 
in Table 131 Although the numbers tend to increase with beam energy, the values presented are 
consistent with each other within the systematic errors. 



Higher quality data would make it possible to derive a more conclusive statement about the 
shape of the curves plotted in Figs. IHl and IHl With the present set of data, usually limited to 
Np above 50, a large part of the centrality curve is missing or smeared by systematic errors. To 
avoid this, we compare Au + Au collisions to p + p (Np=2) at the same energy. 

Figure El shows dNch/dri/{0.5Np) divided by the parameterization plotted in the right panel 
of Fig. 1101 The top panel shows the most central events with Np «350. All points are consistent 
with one demonstrating an agreement of the fit to the data. The points are connected with a line 
for visibility. The middle panel shows results for mid-central events with Np = 100 connected 
with a solid line. The dashed line is the line from the top panel for Np = 350. The points for 
Np = 100 are lower than Np = 350 by a factor of 0.8 — 0.9, over the plotted range of incident 
energies. The lower panel shows p + p data corresponding to Np = 2 measured by several 
experiments. Dashed lines are the same as appear in the upper two panels for Np = 350 and 100 



Table 3. "B/A" ratio and parameter a from the fit to the data. Errors are calculated assuming 
a change of the slope of centrality curves within the limits of the "tilt" errors for PHENIX and 
full errors for the averaged data (Table [2)). 
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Figure 12. The three panels show 
dNch/ dr}/ {.f>Np) divided by the logarithmic 
parameterization from Fig. 1101 The panels 
correspond to Np = 350, 100 and 2 {p+p) from 
top to bottom. Au + Au points are connected 
with lines also shown in lower panels for 
comparison. The Au + Au data is tabulated 
in Table [21 p + p data and parameterizations 
dN/dr] = 2.5 - 0.25/n(siVAr) + Qm'iln{sNNf 
(solid line) and dN/dr] = 0.27ln{sNN) - 0.32 
(dashed line) are taken from |52] . 



and the p + p parameterizations from |51[ 1^ . In the range of RHIC energies these points are 
lower by a factor of 0.65 — 0.75 than the most central. Dotted lines show data from the upper 
two panels. 

These results indicate that the centrality curves normalized to the most central collisions 
have a similar shape for all RHIC energies within the errors of available measurements. 

At lower energies the p+p data show rise and become equal to the most central results. That 
would mean that at the lower SPS and AGS energies the centrality profile is independent of Np. 



At the same time the cross over takes place at AGS energy and it would suggest that below it 
the particle production per participant decreases with Np. There is no measurement supporting 
such a statement. Also, the the information about p + p particle production at low energies is 
recalculated from total number of produced particles or identified particles sets which require 
additional assumptions and systematic uncertainties to be assigned to them. Hence, it is not as 
reliable as at higher energies. 

The data for of the FNAL and BHMC collaborations are shown in Fig. as an indication 
of the trend. 



3.5. Comparison to models 

A variety of models attempting to describe the behavior of Et and Nch as a function of centrality 
at different ^/s^ are available. An updated set of model results were collected from several 
theoretical groups to make a comparison as comprehensive as possible^. 

Figures IT^ through II 51 show the comparison between the existing theoretical models^ and the 
data for 19.6, 130 and 200 GeV. Brief descriptions of the models and their main characteristics 
are given below. 
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Figure 13. dN^h/dr] per pair of participants compared to theoretical models. KLN |58j . 
SSHM EKRT [571, Minijet ^ and JAM |60j. The band shows the range of prediction 
for the Minijet model. 



One of the more commonly used Monte Carlo event generators is HIJING [W\ I55j . This 
model, like several others, uses pQCD for initial minijet production, and the Lund string 
model [HEl for jet fragmentation and hadronization. HIJING also includes jet quenching and 
nuclear shadowing. This type of model typically has two components, a soft part proportional to 
Np and a hard part proportional to Ncoii, which partly motivated the discussion in section 13.4.21 

® For the JAM generator authors extracted data by running the code and assigning Np from the Glauber model 
simulations. 

^ Models are presented as the best fit by the polynomial of the lowest degree which is closer than 1% to any 
theoretical point. The polynomial is plotted in the range where points are provided. 



There are also so-called saturation models which also rely on pQCD and predict that at some 
fixed scale the gluon and quark phase-space density saturates, thus limiting the number of 
produced quarks and gluons. An example of this type of model is EKRT [HZj, which is referred 
to as a final state saturation model. In this paper, comparisons are also made to another parton 
saturation type model, KLN .58j. an initial state saturation model, and also to models related 
to HIJING, namely Minijet [HHl and AMPT |61j . ^MPT is a multiphase transport model, and 
extends HIJING by including explicit interactions between initial minijet partons and also final 
state hadronic interactions. Minijet follows the same two-component model as HIJING but also 
incorporates an energy dependent cut-off scale, similar to the saturation models. JAM jnOj uses 
RQMD and UrQMD inspired ideas for the low energy interactions and above resonance region 
soft string excitation is implemented along the lines of HIJING. 

The other models are listed briefly below. SSHM and SFM did not have a designated short 
identifier, so they were named somewhat arbitrarily here, based on the physics the models 
incorporate. SSHM (Saturation for Semi-Hard Minijet) p7 is also a two-component model: 
pQCD-based for semi-hard partonic interactions, while for the soft particle production it uses 
the wounded nucleon model. DSM j63j . the Dual String Model, is basically the Dual Parton 
Model j64j . with the inclusion of strings. SFM (String Fusion Model) |65j . is a string model 
which includes hard collisions, collectivity in the initial state (string fusion), and rescattering 
of the produced secondaries. Finally, there are the hadronic models, LUCIFER 66). a cascade 
model, with input fixed from lower energy data, and LEXUS [^Tj, a Linear Extrapolation of 
Ultrarelativistic nucleon- nucleon Scattering data to nucleus- nucleus collisions. 

The available model results range from predicting (or postdicting) dNch/dr] at one energy 
to predicting both dNch/dr] and dET/dr] at 19.6, 130 and 200 GeV. The models have varying 
success in reproducing the data. 

Figure El shows that KLN model is one of the most successful at describing the dNch / drj 
centrality dependence for all three energies. However, at y/s^ =19.6 GeV the theoretical curve 
is steeper than the data. This results in a reversed centrality dependence relative to the data for 
the 200 GeV to 19.6 GeV ratio. 5'S'i7M describes the 130 and 200 GeV data weh for centralities 
above Np ~100, which is the approximate limit of applicability for this and other saturation 
models. For the less central events, the models are lower than the data. At 19.6 GeV, the 
model values are significantly higher than the data. The saturation model EKRT describes the 
central points at both energies but overshoots the more peripheral data points and thus does 
not reproduce the general centrality dependence of the data. For the non-saturation models 
included in this figure, Minijet reproduces both the overall scale, as well as the centrality and 
energy dependence of the data rather well. 

The JAM model shown in Fig. and most of the models included in Fig. 1141 provided 
values for all three energies: 19.6, 130 and 200 GeV. JAM is rather successful in describing 
centrality shape at ^/s^ =200 GeV. It also consistent with the measurements at 130 GeV and 
partially at 19.6 GeV. In ratios it underestimates the data. For the Ej- results shown in Fig. 1151 
JAM reproduced the centrality shapes well, but underestimates the y^sj^ dependence of the 
data. SFM is in reasonable agreement with the 130 and 200 GeV data, but gives much larger 
values than the data at 19.6 GeV. ^MPT is overall in good agreement with the data for the two 
higher energies, except for the increasing trend in dNch/dr] at the most peripheral events, which 
is not seen in the data. At the lower energy the N^h centrality behavior is underestimated. 
LEXUS severely overshoots the data for all energies, indicating that nucleus-nucleus effects 
are not correctly accounted for. LUGIFER describes the central points at 130 GeV well, but 
undershoots the less central values at this energy. The HIJING models (version 1.37 and a 
new version with implemented baryon junctions, HIJING B-B) only provide points at 130 and 
200 GeV and are in reasonable agreement with the data at those energies, but generally give 
somewhat lower values. The curves shown include quenching and shadowing HIJING. DSM 



describes 19.6 GeV reasonably well for all centralities, and the more central bins for 130 and 
200 GeV, but overpredicts the values for semi-central and peripheral events. 
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Figure 14. Theoretical models compared to dNch/dr] per pair of participants. SFM |65j . 
AMPT jnU, LEXUS ^,LUCIFER ^,HIJING and DSM (MJ- 



Figure 1131 shows the results for the models that provide data for both dNch/dr] and dEx/dr]. 
For dET/dr], LEXUS and SFM consistently overshoot the data for all energies. In the ratio 
Et/Nch , LEXUS gives values that are too low except at the lowest energy 19.6 GeV. That 
indicates that the hadronization mechanism allows too little energy per particle. The SFM gives 
values that are too large, except for the most peripheral bin, this suggests, that the particles 
are assigned transverse masses that are too large. The HIJING versions and the related AMPT 
model are in reasonable agreement with the data for both dEx/dr] and Ex/Nch^. 

Also shown in Fig. El are the ratios of results at 200 GeV to 19.6 GeV, and 200 GeV to 
130 GeV, for dEx/d-q. These results, especially the comparison of the 200 GeV to 19.6 GeV 
data, is intended to make a more precise check of the ^/s^ dependence of the models. SFM 
fails to describe the 19.6 GeV data and thus can not describe the energy dependence probed by 
these ratios, unlike LEXUS which however does not agree well with the individual data curves 
for 19.6, 130 and 200 GeV. AMPT and the HIJING versions reproduce the values of the ratios 
well, as expected since they are in reasonable agreement with the individual curves. ^MPTand 
HIJING are also successful in describing the Ex/N^h ratio, as illustrated in the lower panels of 
Fig.HSl 

To summarize, most models reproduce at least some of the data fairly well, but most also fail 
in describing all the data. Since the model results typically are given without systematic errors, 
it is not entirely straightforward to quantify the level of agreement or disagreement with the 
data. Qualitatively, the models that are most successful in describing both dEx/drj and dNch/dr] 
in terms of the overall trends, both regarding centrality dependence and energy dependence, are 
AMPT, and the HIJING versions. KLN and Minijet unfortunately do not give information on 

* Note that the HIJING versions available at the time the data were collected and used for predictions were in 
worse agreement with the data This was before energy loss and minijet separation/cut-off scale parameters 
were updated. 
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Figure 15. Theoretical models compared to dEx/dri per pair of participants (upper panels) 
and per produced charged particle (lower panels). SFM [HS], AMPT [^J, LEXUS [HZl and 
HIJING 



dEx/dr] but are successful in describing the dNch/dr] results. The dNch/dr] results thus can 
either be described by the initial state saturation scenario (KLN) or by the mini-jet models that 
need an energy-dependent mini-jet cut-off scale as described in j55, i59j to reproduce the data. 



4. Summary 

This paper presents a systematic study of the energy and centrality dependence of the charged 
particle multiplicity and transverse energy at mid-rapidity at ^/s^ = 19.6, 130 and 200 GeV. 

The yields, divided by the number of participant nucleons, show a consistent centrality 
dependence (increase from peripheral to central) between dEx/dr] and dNch/dij for all energies. 
Furthermore, the increase in the ratio Et/Nch from 19.6 GeV to 200 GeV is consistent with a 
20% increase in {itit) with increasing ^ys^. The ratio Ex/Nch shows only a weak centrality 
dependence at RHIC energies. 

For the ^/s^ dependence, comparisons were made not only among RHIC results but 
also including data from lower energy fixed-target experiments at SPS, AGS and SIS. A 
phenomenological fit, scaling logarithmically with y^s^, describes both dEx/dr] and dNch/dr] 
well for all energies and for the most central collisions. 



Using the fit results, one can separate two regions with different particle production 
mechanisms. The region below SPS energy is characterized by a steep increase in Ex/Nch ^{rriT) 
with ^/s^, whereas for the energies above SPS Et/Nch is found to be weakly dependent on 

Within the systematic errors of the measurements the shape of the centrality curves of 
dNch/dr]/{0.5Np) vs. Np were found to be the same in the range of RHIC energies and to 
scale with ln(y^ijV7v)- The same trend must be true for Et because Et/Nch has a very weak 
centrality dependence. 

Based on the dET/dr] measurements, the Bjorken energy density estimates were performed 
and esj ■ t was determined to be 5.4 it 0.6 GeV/(/m^ • c) at ^/s^ = 200 GeV for the most 
central bin. This is in excess of what is believed to be sufficient for a phase transition to the 
new state of matter. The energy density increases by about a factor of 2.6 from SPS to RHIC 
energy. 

Finally, a comparison between the RHIC dNch/dr] and dET/dr] data and a collection of 
models was performed. A few models, notably HIJING and AMPT, reproduce both dET/dr] 
and dNch/dr] rather well for several energies. 
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